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\ ABSTRACT 

' . We present a cosmic shear analysis of the 100 square degree weak lensing survey, combining 

p~ ' data from the CFHTLS-Wide, RCS, VIRMOS-DESCART and GaBoDS surveys. Spanning 

| ~ 100 square degrees, with a median source redshift z ~ 0.78, this combined survey allows 

— . i us to place tight joint constraints on the matter density parameter fi m , and the amplitude of 

'"^ \ the matter power spectrum erg, finding crs(^f|) ' 59 = 0.84 ± 0.05. Tables of the measured 

shear correlation function and the calculated covariance matrix for each survey are included 
O ' as Supplementary Material to the online version of this article. 

\ The accuracy of our results are a marked improvement on previous work owing to three 

important differences in our analysis; we correctly account for sample variance errors by in- 
cluding a non-Gaussian contribution estimated from numerical simulations; we correct the 
measured shear for a calibration bias as estimated from simulated data; we model the red- 
shift distribution, n(z), of each survey from the largest deep photometric redshift catalogue 
currently available from the CFHTLS-Deep. This catalogue is randomly sampled to repro- 
duce the magnitude distribution of each survey with the resulting survey dependent n(z) 
parametrised using two different models. While our results are consistent for the n(z) models 
tested, we find that our cosmological parameter constraints depend weakly (at the 5% level) 
on the inclusion or exclusion of galaxies with low confidence photometric redshift estimates 
(z > 1.5). These high redshift galaxies are relatively few in number but contribute a sig- 
nificant weak lensing signal. It will therefore be important for future weak lensing surveys to 
obtain near-infra-red data to reliably determine the number of high redshift galaxies in cosmic 
shear analyses. 

Key words: cosmology: cosmological parameters - gravitational lenses - large-scale structure 
of Universe - observations 



1 INTRODUCTION 

Weak gravitational lensing of distant galaxies by intervening matter 
provides a unique and unbiased tool to study the matter distribution 
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of the Universe. Large scale structure in the Universe induces a 
weak lensing signal, known as cosmic shear. This signal provides 
a direct measure of the projected matter power spectrum over a 
redshift range determined by the lensed sources and over scales 
ranging from the linear to non-linear regime. The unambiguous in- 
terpretation of the weak lensing signal makes it a powerful tool for 
measuring cosmological parameters which compl ement those from 
other probes such a s, CMB anisotropi es JSpergeletaii r2006) and 
type la supernovae jAstieretalll2006h . 

Cosmic shear has only recently become a viable tool for 
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observational cosmology, with the fi r st measurements reported 
simultaneously in 2000 dBacon etail |2000|: iKaiser et all hood : 



IVan Waerbeke et all l200d Iwittman et all l200Ch . The data these 



early studies utilised were not optimally suited for the extraction 
of a weak lensing signal. Often having a poor trade off between 
sky coverage and depth, and lacking photometry in more than one 
colour, the ability of these early surveys to constrain cosmology 
via weak lensing was limited, but impressive based on the data 
at hand. Since these early results several large dedicated surveys 



petitive constraints 


on cosmological parameters ( Hoekstra et al 


2002al: IBacon et all 


20031: IJarvis et all |2003|: iBrown et all 


2003; 


Hamana et al., 2003 


; Massey et al., 200a iRhodes et al., 


2004 


Van Waerbeke et al., 


2005|;[Hevmans et al., 2005; Semboloni et al. 


20061; iHoekstra et al 


. 2006; Schrabback et alll2006t). With the re- 



cent first m easurements of a changing lensing signal as a function 
of redshift (IBacon et all 120051 ; Iwittmanl 120051; ISemboloni et all 
2006; Mas s ev et all l2007bl), and the fir s t weak lensing con- 



straints on dark energy djarvis et all 120061 ; iKitching et all |2006| ; 
IHoekstra etlull2006l : lschimd et al.Ll2007h . the future of lensing is 



promising. 

There are currently several 'next generation' surveys being 
conducted. Their goal is to provide multi-colour data and excellent 
image quality, over a wide field of view. Such data will enable weak 
lensing to place tighter constraints on cosmology, breaking the de- 
generacy between the matter density parameter Q m and the ampli- 
tude of the matter power spectrum <xs and allowing for competitive 
constraints on dark energy. In this paper, we present an analysis of 
the 100 deg 2 weak lensing survey that combines data from four of 
the largest surveys analysed to date, including the wide component 
of the Canada-France-Hawa ii Telescope Legacy Survey (CFHTLS- 
Wide, Hoekstra et al. , 2006), the Garching-Bonn Deep Survey (Ga- 
BoDs, Hetterscheidt et al., 2006), the Red-Sequence Cluster Sur- 
vey (RCS. IHoekstra et aul2002al). and the V IRMOS-DESCART 
survey (VIRMOS. IVan Waerbeke et alll2005l) . These surveys have 
a combined sky coverage of 113 deg 2 (96.5 deg 2 , after masking), 
making this study the largest of its kind. 

Our goal in this paper is to provide the best estimates of cos- 
mology currently attainable by weak lensing, through a homoge- 
neous analysis of the major data sets available. Our analysis is dis- 
tinguished from previous work in three important ways. For the first 
time we account for the effects of non-Gaussian contributions to the 
analytic estimate of the sample varianc e - sometimes referre d to 
as cosmic v ariance - covariance ma trix dSchneider et alll2002l) . as 
described in ISemboloni etai] d2007l). Results from the Shear TEst- 
ing Pr ogramme (STEP, see iMassev et all l2007al ; iHevmans et all 
2006a) are used to correct for calibration error in our shear mea- 
surement methods, a marginalisation over the uncertainty in this 
correction is performed. We use th e largest deep pho tometric red- 
shift catalogue currently available dllbert et all 12006'). which pro- 
vides redshifts for ~ 500000 galaxies in the CFHTLS-Deep. Addi- 
tionally, we properly account for sample variance in our calculation 
of the redshift distribu t ion, an important source of error discussed 
in lVan Waerbeke et alj d2006h . 

This paper is organised as follows: in §2 we give a short 
overview of cosmic shear theory, and outline the relevant statis- 
tics used in this work. We describe briefly each of the surveys used 
in this study in §3. We present the measured shear signal in §4. 
In §5 we present the derived redshift distributions for each survey. 
We present the results of the combined parameter estimation in §6, 
closing thoughts and future prospects are discussed in §7. 



2 COSMIC SHEAR THEORY 

We briefly describe here the notations and statistics used in our cos- 
mic shear analysis. For detailed reviews o f weak lensing theory the 
reader is re ferred tolMunshi et alJd2~006t) . lBartelmann & Schneideil 
(2001), and lSchneider et al.ldl998T) . The notations used in the latter 
are adopted here. The power spectrum of the projected density field 
(convergence ft) can be written as 



4c 4 



dw 

a 2 (w) 



dw n(w ) 



f K (w) 
2 



(1) 



where Ho is the Hubble constant, /k{w) is the comoving angu- 
lar diameter distance out to a distance w (wh is the comoving 
horizon distance), a(w) is the scale factor, and n[ui(z)] is the red- 
shift distribution of the sources (see §5). P3D is the 3-dimensional 
mass power spectrum computed from a n on-linear estimation of the 
dark matter cluste ring (see for example IPeacock & Doddsl 1 19961 : 
ISmith et all l2003h . and k is the 2-dimensional wave vector per- 
pendicular to the line-of-sight. P30 evolves with time, hence it's 
dependence on the co-moving radial coordinate w. 

In this paper we focus on the shear correlation function statis- 
tic £. For a galaxy pair separation 9, we define 



= {7t(rht(r + 6)) + <rtr(ry/r(r + 9)). 
= (7t(r)-ft(r + 0))-<7r(r)7r(r + 0)>, 



(2) 



where the shear 7 = (7t,7r) is rotated into the local frame of 
the line joining the centres of each galaxy pair separated by 9. The 
tangential shear is 74, while j r is the rotated shear (having an angle 
of -| to the tangential component). The shear correlation function 
£+ is related to the convergence power spectrum through 



1 

2^ 



dk kP K (k)J (k9), 



(3) 



where Jo is the zeroth order Bessel function of the first kind. A 
quantitative measurement of the lensing amplitude and the sys- 
tematics is obtained by splitting the signal into its curl-free (E- 
mode) and curl (B-mode) components respectively. This method 
has been advocated to help the measurement of t he intrinsic align- 
ment contam ination in the weak lensing signal dCrittenden et all 
I200lll2002l) , but it is also an effici ent measure of the residual sys- 
tematics from the PSF correction dPen etall 120021) . 

The E and B modes derived from the shape of galaxies are 
unambiguously defined only for the so-called aperture mass vari- 
ance (Mf p ), which is a weighted shear variance within a cell 
of radius 9 C . The cell itself is defined as a compensated filter 
dSchneider et all 1998), such that a constant convergence ft gives 
Af ap = 0. (Mf p ) can be rewritten as a function of the tangential 
shear 74 if we express 7 = (7t,7r) in the local frame of the line 
connecting the aperture centre to the galaxy. (Map) is given by: 



288 



dk 

^P K (k)[J 4 (k9 c )} 2 . 



(4) 



The B-mode (M ap )x is obtained by replacing 74 with 7 r . The 
aperture mass is insensitive to the mass sheet degeneracy, and there- 
fore it provides an unambiguous splitting of the E and B modes. 
The drawback is that aperture mass is a much better estimate of 
the small scale power than the large scale power which can be seen 
from the function Ji(k9 c ) in Eq.«(4j which peaks at k9 c ~ 5. Es- 
sentially, all scales larger than a fifth of the largest survey scale re- 
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main inaccessible to M ap . The large-scale part of the lensing signal 
is lost by M a p, while the remaining small-scale fraction is difficult 
to interpret beca use the strongly non- l inear p ower is difficult to pre- 
dict accurately JVan Waerbeke et all 120021) . It is therefore prefer- 
able to decompose the shear correlation function into its E and B 
modes, as it is a much deeper probe of the lin ear regime. 
Following lcrittenden et af] j200 11 12002). we define 



desig ned to match the SDSS photometric system jFukugita et all 



L'UD = ( u))+ l / ^£_(tf)-120 2 



The E and B shear correlation functions are then given by 



(6) 



2 ' 2 

In the absence of systematics, £b — and £e = £+ (Eq|3j. In 
contrast with the M ap statistics, the separation of the two modes 
depends on the signal integrated out with the scales probed by all 
surveys. One option is to calculate Eq.l(6j using a fiducial cosmol- 
ogy to compute £_ on scales 8 —-> oo. As shown in lHevmans et al.l 
(2005), chan ges in the choice of fiducial cosmology do not sig- 
nificantly affect the results, allowing this statistic to be used as a 
diagnostic tool for the presence of systematics. This option does 
however prevent us from using £e to constrain cosmology. As the 
statistical noise on the measured £e is ~ \/~2 smaller than the sta- 
tistical noise on the measured £ + , there is a more preferable al- 
ternative. If the survey size is sufficiently large in comparison to 
the scales probed by ^b, we can consider the unknown integral to 
be a constant that can be calibrated using the aperture mass B- 
mode statistic (M ap (A0))j_. A range A9 of angular scales where 
(Affp (A0))j_ ~ ensures that the J3-mode of the shear correla- 
tion function is zero as well (within the error bars), at angular scales 
~ In this analysis we have calibrated £,e,b for our survey data 
in this manner. This alternative method has also been verified by re- 
calculating £e,b using our final best-fit cosmology to extrapolate 
the signal. We find the two methods to be in very close agreement. 



3 DATA DESCRIPTION 

In this section we summarise the four weak lensing surveys 
that fo rm the 100 d eg 2 w eak lensing survey: CFHTLS- 
Wide dHoekstra etall l200d), Ga BoDs dHetterscheidt et all 



l2006t). RCS dHoekstra et all 12002 ah . and VIRMOS-DESCART 
dVan Waerbeke et all 120051) , see Table [T] for an overview. Note 
that we have chosen not to include the CFHTLS-Deep data 
dSemboloni etail l2006h in this analysis since the effective area 
of the analysed data is only 2.3 deg 2 . Given the breadth of the 
parameter constraints obtained from this preliminary data set there 
is little value to be gained by combining it with four surveys which 
are among the largest available. 

3.1 CFHTLS-Wide 

The Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) 
is a joint Canadian-French program designed to take advantage 
of Megaprime, the CFHT wide-field imager. This 36 CCD mo- 
saic camera has a 1 x 1 degree field of view and a pixel scale of 
0.187 arcseconds per pixel. The CFHTLS has been allocated ~ 450 
nights over a 5 year period (starting in 2003), the goal is to com- 
plete three independent survey components; deep, wide, and very 
wide. The wide component will be imaged with five broad-band 
filters: u* , g\ r', i , z' ■ With the exception of u* these filters are 



1996). The three wide fiel ds will span a total of 170 deg 2 on com- 
pletion. In this study, as in lHoekstra et al. (2006), we use data from 
the Wl and W3 fields that total 22 deg 2 after masking, and reach 
a depth of i' = 24.5. For detailed information con cerning data re- 
ductio n and object analysis the reader is referred to lHoekstra et al.l 
d2006l) . and lVan Waerbeke et alj J2OO2I) . 



3.2 GaBoDS 

The Garching-Bonn Deep Survey (GaBoDS) data set was obtained 
with the Wide Field Imager (WFI) of the MPG/ESO 2.2m tele- 
scope on La Silla, Chile. The camera consists of 8 CCDs with a 34 
by 33 arcminute field of view and a pixel scale of 0.238 arcseconds 
per pixel. The total data set consists of 52 statistically independent 
fields, with a total effective area (after trimming) of 13 deg 2 . The 
magnitude limits of each field ranges between J?vega = 25.0 and 
Evega = 26.5 (this corresponds to a 5<r detection in a 2 arcsec- 
ond apertu re). To obtain a ho mogeneous depth across the survey 
we follow iHetterscheidt et all ( 120061) performing our lensing anal- 
ysis with only those objects which lie in the complete magnitude 
interval R £ [21.5, 24.5]. For further details of the data set, and 
technical i nformation regarding the r eduction pipeline the reader is 
referred to lHetterscheidt et all ( 120061) . 



3.3 VIRMOS-DESCART 

The DESCART weak lensing project is a theoretical and observa- 
tional program for cosmological weak lensing investigations. The 
cosmic shear survey carried out by the DESCART team uses the 
CFH12k data jointly with the VIRMOS survey to produce a large 
homogeneous photometric sample in the BVRI broad band filters. 
CFH12k is a 12 CCD mosaic camera mounted at the Canada- 
France-Hawaii Telescope prime focus, with a 42 by 28 arcminute 
field of view and a pixel scale of 0.206 arcseconds per pixel. 
The VIRMOS-DESCART data consist of four uncorrelated patches 
of 2 x 2 deg 2 separated b y more than 40 degrees. Here, as in 
IVan Waerbeke et al.l (2005), we use Jab -band data from all four 
fields, which have an effective area of 8.5 deg 2 after masking, and a 
limiting magnitude of Jab = 24.5 (this corresponds to a 5a detec- 
tion in a 3 arcsecond aperture). Tech nica l details of the data set ar e 
given in lVan Waerbeke et al.l d200ll) and iMcCracken et al.l d2003h . 
while an overview of the VIRMOS-DESCART survey is given in 
iLe Fevreetai] ( l2004l) . 



3.4 RCS 

The Red-Sequence Cluster Survey (RCS) is a galaxy cluster survey 
designed to provide a large sample of optically selected clusters 
of galaxies with redshifts between 0.3 and 1.4. The final survey 
covers 90 deg 2 in both Rc and z , spread over 22 widely sepa- 
rated patches of ~ 2.1 x 2.3 degrees. The northern half of the 
survey was observed using the CFH12k camera on the CFHT, and 
the southern half using the Mosaic II camera mounted at the Cerro 
Tololo Inter-American Observatory (CTIO), Victor M. Blanco 4m 
telescope prime focus. This camera has an 8 CCD wide field imager 
with a 36 by 36 arcminute field of view and a pix el scale of 0.27 
arcseco nds per pixel. The data studied here, as in iHoekstra et al.l 
d2002al) . consist of a total effective area of 53 deg 2 of masked imag- 
ing data, spread over 13 patches, with a limiting magnitude of 25.2 
(corresponding to a 5<r point source depth) in the Rc band. 43 deg 2 
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Table 1. Summary of th e published results from each survey used in this study. The values for <rg correspond to f2 m = 0.24, and are given for the 
IPeacock & DoddsMl99fj|) method of calculating the non-linear evolution of the matter power spectrum. The statistics listed were used in the previous analyse s; 
the shear correlation functions (,e,b (Eq.|6), the aperture mass statistic (M 2 p ) (Eq.[4) and the top-hat shear variance g ) (see lVan Waerbeke et alj feoOS)). 



CFHTLS-Wide GaBoDS RCS VIRMOS-DESCART 

Area(deg 2 ) 22.0 13.0 53 8.5 

N flcUs 2 52 13 4 

Magnitude Range 21.5 < i' < 24.5 21.5 < R < 24.5 22 < R c < 24 21 < 7 AB < 24.5 

< Zsource > 0.81 0.78 0.6 0.92 

^median 0.71 0.58 0.58 0.84 

Previous Analysis Hoekstra et al. (2006) Hetterscheidt et al. (2006) Hoekstra et al. (2002a) Van Waerbe ke et al. (2005) 

cr 8 (n m = 0.24) 0.99 ±0.07 0.92 ±0.13 0.98 ±0.16 0.96 ±0.08 

Statistic £ B>B (0) (M 2 ap )(e) {Ml p ){9) <7|, B >(0) 



were taken with the CFHT, the remaining 10 deg 2 were taken with 
the Blanco telescope. A detailed description of the data reduction 
and object analysis is described in Hoekst ra et al. (2002b), to which 
we refer for technical details. 



the error on the calibration bias, which expresses the range of ad- 
missible scalings of Tables containing the measured correlation 
functions and their corresponding covariance matrices are included 
as Supplementary Material to the online version of this article. 



4 COSMIC SHEAR SIGNAL 



Lensing shear has been measured for each sur vey using the com- 
mon KSB galaxy shape measurement method dKaiser et all 1 19951 ; 
Luppino & Kaiseil [l997l ; iHoekstra et all 1 19981) . The practical ap- 
plication of this method to the four surveys di ffers slightly, a s tested 
by the Shear TEsting Programme; STEP ( Ma ssev et all [2007a; 
Heymans et al., 2006a), r esulting in a small calib ration bias in the 
measurement of the shear. lHevmans et al (2006a) define both a cal- 
ibration bias (m) and an offset of the shear (c), 



In 



7true = Tl7truo + C, 



(7) 



where 7 mca s is the measured shear, 7t ruc is the true shear signal, 
and m and c are determined from an analysis of simulated data. 
The value of c is highly dependent on the strength of the PSF dis- 
tribution, therefore the value determined through comparisons with 
simulated data can not be easily applied to real data whose PSF 
strength varies across the image. Since an offset in the shear will 
appear as residual B-modes, we take c = 0, and focus on the more 
significant and otherwise unaccounted for shear bias (m). Eq.Q 
then simplifies to 



7tr 



= 7mcas("J+l) • 



(8) 



The calibration bias for the C FHTLS-Wide RCS an d VIRMOS- 
DESCART (HH analysis in iMassev et all J2007al) ) is m = 
-0.0017 ± 0.0088, and ro = 0.038 ± 0.026 for GaBoDS (MH 
analysis in lMassev et alj J2007ah ). 

We measure the shear correlation functions £e(8) and (,b(9) 
(Eq.(6]( for each survey, shown in Figure[T] As discussed in §|2] these 
statistics are calibrated using the measured aperture mass B-mode 
{Mlp) X • T he l(j errors on £e incl ude statistical noise, computed as 
described in Schneider et al. (2002), non-Gaussian sample variance 
(see §6 for details), and a systematic error added in quadrature. The 
systematic error for each angular scale 6 is given by the magnitude 
of £,b{0). The 1<t errors on £s are statistical only. We correct for 
the calibration bias by scaling £e by (1 + m)~ 2 , the power of neg- 
ative 2 arising from the fact that £s is a second order shear statistic. 
When estimating parameter constraints in §6.2 we marginalise over 



5 REDSHIFT DISTRIBUTION 

The measured weak lensing signal depends on the redshift distribu- 
tion of the sources, as seen from Eq.Q}. Past weak l ensing studies 
(e.g.jHoekstra et al.1 j2006h;ISemboloiii et alj j2006l) ; lMassev et all 
j2005h ; IVan Waerbeke et al.1 j2005h ) have used the Hubble Deep 
Field (HDF) photometric redshifts to estimate the shape of the red- 
shift distribution. Spanning only 5.3 arcmin 2 , the HDF suffe rs from 
sample variance as described in lVan Waerbeke et alj ( 1200a) , where 
it is also suggested that the HDF fields may be subject to a selec- 
tion bias. This sample variance of the measured redshift distribu- 
tion adds an additional error to the weak lensing analysis that is not 
typically accounted for. 

In this study we use t he largest d e ep ph otometric redshift cat- 
alogue in existence, from lllbert et al.1 J2006l) . who have estimated 
redshifts on the four Deep fields of the T0003 CFHTLS release. 
The redshift catalogue is publicly available at terapix.iap.fr . The 
full photometric catalogue contains 522286 objects, covering an 
effective area of 3.2 deg 2 . A set of 3241 spectroscopic redshifts 
with ^ z «C 5 from the VIRMOS VLT Deep Survey (VVDS) 
were used as a calibration and training set for the photometric 
redshifts. The resulting photometric redshifts have an accuracy of 
o-( Zphot -z apcc )/(i+ Zapoc ) = 0.043 for i' ab = 22.5 - 24 , with a frac- 
tion of catastrophic errors of 5.4%. lllbertetal.l (2006) demonstrate 
that their derived redshifts work best in the range 0.2 ^ z ^ 1.5, 
having a fraction of catastrophic errors of ~ 5% in this range. The 
fraction of catastrophic errors increases dramatically at z < 0.2 
and 1.5 < z < 3 reaching ~ 40% and ~ 70% respectively. This 
is explained by a degeneracy between galaxies at z spcc < 0.4 and 
1.5 < Zphot < 3 due to a mismatch between the Balmer break and 
the i ntergalactic Lyman- alpha forest depression. 

IVan Waerbeke et al. ( 2006) estimate the expected sampling er- 
ror on the average redshift for such a photometric redshift sample 
to be ~ 3%/Vi = 1-5%, where the factor of comes from the 
four independent CFHTLS-Deep fields, a great improvement over 
the ~ 10% error expected for the HDF sample. 
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Figure 1. E and B modes of the shear correlation function £ (filled and open points, respectively) as measured for each survey. Note that the la errors on the 
B-modes include statistical noise, non-Gaussian sample variance (see §6) and a systematic error given by the magnitude of the ZJ-mode. The la error on the 
B-modes is statistical only. The results are presented on a log-log scale, despite the existence of negative B-modes. We have therefore collapsed the infinite 
space between 10 — 7 and zero, and plotted negative values on a separate log scale mirrored on 10~ 7 . Hence all values on the lower portion of the graph are 
negative, their absolute value is given by the scaling of the graph. Note that this choice of scaling exaggerates any discrepancies. The solid lines show the best 
fit ACDM model for Q m = 0.24, h = 0.72, T = hD. m , erg given in Table l4l and n(z) given in Table |2] The latter t wo being chosen fo r the case of the high 
confidence redshift calibration sample, an rt(z) modeled by Eq. jlOh and the non-linear power spectrum estimated by Smith et al. (20031). 



Table 2. The best fit model parameters for the redshift distribution, given by either Eq.fiOl or Eq.{TT). Models are fit using the full calibration sample of 
lllbert et al.l 12006) 0.0 $C Zp ^ 4.0, and using only the high confidence region 0.2 ^ z v ^ 1.5. i s tne reduced % 2 statistic, (z), and z m are the average 
and median redshift respectively, calculated from the model on the range 0.0 ^ z ^ 4.0. 
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5.67 


0.819 


4.418 


0.800 


1.201 


0.939 


0.808 


1.84 




VIRMOS-DESCART 


1.045 


1.445 


0.767 


0.909 


0.816 


9.95 


0.703 


5.000 


1.763 


2.042 


0.960 


0.864 


2.40 
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5.1 Magnitude Conversions 

In order to estimate the redshift distributions of the surveys, we cal- 
ibrate the magnitude distribution of the photometric redshift sample 
(hence forth the z p sample) to that of a given survey by converting 
the CFHTLS filter set u'g'r'i'z to 7 A b magnitudes for VIRMOS- 
DESCART and .Rvega magnitudes for RCS and GaBoDS. We 
e mploy the linear rela t ionshi ps between different filter bands given 
bv lBlanton & Roweisl J2007h ; 

Jab = i' - 0.0647 - 0.7177 [{% - z) - 0.2083] , (9) 

flvEGA = r - 0.0576 - 0.3718 [(r - i') - 0.2589] - 0.21. 

These conversions were estimated by fitting spectral energy distri- 
bution templates to data from the Sloan Digital Sky Survey. They 
are considered to be accurate to 0.05 mag or better, resulting in an 
error on our estimated median redshifts of at most 1%, which is 
small compared to the total error budget on the estimated cosmo- 
logical parameters. 

5.2 Modeling the redshift distribution n(z) 

The galaxy weights used in each survey's lensing analysis result 
in a weighted source redshift distribution. To estimate the effec- 
tive n(z) of each survey we draw galaxies at rand om from the z v 
sampl e, using the method described in §6.5.1 of IWall & Jenkins! 
( 2003) to reproduce the shape of the weighted magnitude distribu- 
tion of each survey. A Monte Carlo (bootstrap) approach is taken 
to account for the errors in the photometric redshifts, as well as the 
statistical variations expected from drawing a random sample of 
galaxies from the z p sample to create the redshift distribution. The 
process of randomly selecting galaxies from the z p sample, and 
therefore estimating the redshift distribution of the survey, is re- 
peated 1000 times. Each redshift that is selected is drawn from its 
probability distribution defined b y the ±1 and ±3a errors (since 
these are the error bounds given in Ilb ert et alj j2006T) 's catalogue). 
The sampling of each redshift is done such that a uniformly ran- 
dom value within the la error is selected 68% of the time, and a 
uniformly random value within the three sigma error (but exterior 
to the la error) is selected 32% of the time. The redshift distribu- 
tion is then defined as the average of the 1000 constructions, and 
an average covariance matrix of the redshift bins is calculated. 

At this point sample variance and Poisson noise are added 
to the diagonal elements of the covariance matrix, following 
Van Waerbeke et al. (2006). They provide a scaling relation be- 
tween cosmic (sample) variance noise and Poisson noise where 
Csampie/cpoisson is given as a function of redshift for different 
sized calibration samples. We take the curve for a 1 deg 2 survey, 
the size of a sing le CFHTLS-Deep field, and divide by y/l since 
there are four independent 1 deg 2 patches of sky in the photometric 

redshift sample. 

Based on the photometric redshift sample of lllbert et all 

j2006l) we consider two redshift ranges, the full range of redshifts 
in the photometric catalogues 0.0 ^ z ^ 4.0, and the high confi- 
dence range 0.2 ^ z ^ 1.5. The goal is to assess to what extent 
this will affect the redshift distribution, and -in turn- the derived 
parameter constraints. 

The shape of the normalised redshift distribution is often as- 
sumed to take the following form: 




where a, /3, and zo are free parameters. If the full range of pho- 
tometric redshifts is used the shape of the redshift distribution is 
poorly fit by Eq.dlOt, this is a result of the function's exponential 
drop off which can not accommodate the number of high redshift 
galaxies in the tail of the distribution (see Fig.l[2j). We therefore 
adopt a new function in an attempt to better fit the normalised red- 
shift distribution, 

n(z)=N-^-, (11) 
where a, b, and c are free parameters, and TV is a normalising factor, 

N ={J dz '—) • (12) 

The best fit model is determined by minimising the gener- 
alised chi square statistic: 

X = (di - m»)C _1 (di - m;) T , (13) 

where C is the covariance matrix of the binned redshift distribution 
as determined from the 1000 Monte Carlo constructions, di is the 
number count of galaxies in the i th bin from the average of the 1000 
distributions, and nn is the number count of galaxies at the center 
of the i th bin for a given model distribution. For each survey we 
determine the best fit for both Eq. d lQt and Eq.d lib. In either case 
we consider both the full range of redshifts and the high confidence 
range, the results are presented in Table|2] 

Figure [2] shows the best fit models to the CFHTLS-Wide sur- 
vey, when the full range of photometric redshifts is considered (left 
panel), Eq.il It clearly fits the distribution better than Eq.lllOt, hav- 
ing reduced \ 2 statistics of 2.41 and 8.94 respectively. When we 
consider only the high confidence redshifts (right panel of Figure(2]l 
both functions fit equally well having reduced \ 2 statistics of 1.65 
and 1 .63 for Eq.J 1 It and Eq. dlOt respectively. Note that the redshift 
distribution as determined from the data (histograms in Figure|2j is 
non-zero at z > 1.5 for the high confidence photometric redshifts 
(0.2 ^ z p 1.5) because of the Monte Carlo sampling of the 
redshifts from their probability distributions. 



6 PARAMETER ESTIMATION 
6.1 Maximum likelihood method 

We investigate a six dimensional parameter space consisting of the 
mean matter density f2 m , the normalisation of the matter power 
spectrum as, the Hubble parameter h, and n(z) the redshift distri- 
bution parametrised by either a, (3, and zo (EqllOt or a, b, and c 
(Eqll It. A flat cosmology (fi m + Qa = 1) is assumed throughout, 
and the shape parameter is given by V — fi m h. The default priors 
are taken to be fi m G [0.1, 1], a 8 G [0.5, 1.2], and h G [0.64, 0.8] 
with the latter in agreem ent with the findings of the HST key project 
dFreedman et all uOOl). The priors on the redshift distribution were 
arrived at using a Monte Carlo technique. This is necessary since 
the three parameters of either Eq.lllOt or Eq.J 1 It are very degen- 
erate, hence simply finding the 2 or 3 a levels of one parameter 
while keeping the other two fixed at their best fit values does not 
fairly represent the probability distribution of the redshift param- 
eters. The method used ensures a sampling of parameter triplets 
whose number count follow the 3-D probability distribution; that 
is 68% lie within the la volume, 97% within the 2a volume, etc. 
We find 100 such parameter trios and use them as the prior on the 
redshift distribution, therefore this is a Gaussian prior on n(z). 
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Figure 2. Normalised redshift distribution for the CFHTLS-Wide survey, given by the histogram, where the error bars include Poisson noise and sample 
variance of the photometric redshift sample. The dashed curve shows the best fit for Eq. tlOt . and the solid curve for Eq.jTT}. Left: Distribution obtained if 
all photometric redshifts are used 0.0 ^ z $C 4.0, and \ 2 is calculated between 0.0 $C z $C 2.5. Right: Distribution obtained if only the high confidence 
redshifts are used 0.2 $C z $C 1.5, and \ 2 is calculated on this range. The existence of counts for z > 1.5 is a result of drawing the redshifts from their full 
probability distributions. 



Given the data vector which is the shear correlation func- 
tion (£e of Eq[6]l as a function of scale, and the model prediction 
m(Sl m , <7g, h, n(z)) the likelihood function of the data is given by 



exp 



-(|-m)C- 1 (^-m) 



(14) 



where n is the number of angular scale bins and C is the n x n 
covariance matrix. The shear covariance matrix can be expressed 
as 



(15) 



where fn is the mean of the shear £i at scale i, and angular brackets 
denote the average over many independent patches of sky. To obtain 
a reasonable estimate of the covariance matrix for a given set of 
data one needs many independent fields, this is not the case for 
either the CFHTLS-Wide or VIRMOS-DESCART surveys. 

We opt to take a consistent approach for all 4 data sets by de- 
composing the shear covariance matrix as C = C n + Cb + C s , 
where C n is the statistical noise, Cb is the absolute value of the 
residual B-mode and C s is the sample variance covariance matrix. 
C n can be measured directly from the data, it represents the statisti- 
cal noise inherent in a finite data set. Cb is diagonal and represents 
the addition of the B-mode in quadrature to the uncertainty, this 
provides a conservative limit to how well the lensing signal can be 
determined. 

For Gaussian sample variance the m atrix C s can be com- 
puted according to ISchneider et~ai1 ( l2002h . assuming an effective 
survey area A c g, an effective number density of galaxies n e g, 
and an intrinsic ellipt icity dispersion a e (see Table [3). However, 
Schnei der et al. assume Gaussian statistics for the fourth or- 

der moment of the shear correlation function -a necessary simpli- 
fication to achieve an analytic form for the sample variance covari- 
ance m atrix. We use the calibration presented by Sembolo ni et al.l 
d2007l) . estimated from ray tracing simulations, to account for non- 
Gaussianities. Their work focuses on the following quantity: 



•F(0i,02) = 



cr s K+;^i,^) 



(16) 



Table 3. Data used for each survey to calculate the Gaussian contribution to 
the sample variance covariance matrix. A e ff is the effective area, n c ff the 
effective galaxy number density, and cr E the intrinsic ellipticity dispersion. 





(deg 2 ) 


"off 

(arcmin~ 2 ) 




CFHTLS-Wide 


22 


12 


0.47 


RCS 


53 


g 


0.44 


VIRMOS-DESCART 


8.5 


15 


0.44 


GaBoDS 


13 


12.5 


0.50 



where $2) is the ratio of the sample variance covari- 

ance measured from N-body simulations (C™ cas (£+; $1, $2)) to 
that expected from Gaussian effects alone (C? aus (£+; $1, $2))- 
It is found that J-('&i,'&2) increases significantly above unity at 
scales smaller than ~ 10 arcminutes, increasing with decreasing 
scale, it reaches an order of magnitude by ~ 2 arcminutes. The 
parametrised fit as a function of mean source redshift is given by, 



^(#1,^2) 
Pi(z) 



Pi(z) 



[i9 2 tf 2 ] P2(z) 



(17) 



16.9 

^0.95 " 

1.62 z~ 



2.19, 



exp(- 



0.03. 



A fiducial model is required to calculate the Gaussian co- 
variance matrix, it is taken as S7 m = 0.3, A = 0.7, as = 0.8, 
h — 0.72, and the best fit n(z) model (see Table[2]l. We then use 
the above prescription to account for non-Gaussianities, increasing 
C s everywhere T is above unity. For each survey the total covari- 
ance matrix (C = C n + C s + Cb) is included as Supplementary 
Matrial to the online version of this article. 

To test this method we compare our analytic covariance ma- 
trix for GaBoDS with tha t found by measuring it from the data 
dHetterscheidt et alT . l2006t) . Since GaBoDS images 52 independent 
fields it is possible to obtain an estimate of C directly from the data. 
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Figure 3. Joint constraints on eg and f2 m from the 100 deg weak lensing 
survey assuming a flat ACDM cosm ology and adopting the non-linear mat- 
ter power spectrum of Smith et alj fe003l) . The redshift distributio n is esti- 
mated from the high confidence photometric redshift catalogue l llb ert et all 
l2006h . and modeled with the standard functional form given by Eq.flO). 
The contours depict the 0.68, 0.95, and 0.99% confidence levels. The mod- 
els are marginalised, over h = 0.72 ± 0.08, shear calibration bias (see §4) 
with uniform priors, and the redshift distribution with Gaussian priors (see 
§6.1). Similar results are found for all other cases, as listed in Tabled 
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Figure 4. Joint constraints on <rg and n m for each survey assuming a flat 
ACDM cosmology and adopting the non-linear matter power spectrum of 
ISmifhetalJl2003l) . The redshift distribution is estimated from the high con- 
fidence photometric redshift catalogue Hlbert et all |2006J), and modeled 
with the standard functional form given by Eq. dlOt . The contours depict 
the 0.68, 0.95, and 0.99% confidence levels. The models are marginalised, 
over h = 0.72 ± 0.08, shear calibration bias (see §4) with uniform priors, 
and the redshift distribution with Gaussian priors (see §6.1). Similar results 
are found for all other cases, as listed in Table|4] 



The contribution from the B-modes (Cs) is not added to our ana- 
lytic estimate in this case, since its inclusion is meant as a conser- 
vative estimate of the systematic errors. We find a median percent 
difference along the diagonal of ~ 15% , which agrees well wi th 
the accuracy obtained for simulated data JSemboloni et alll2007h . 



6.2 Joint constraints on erg and fi m 

For the combined survey we place joint constraints on erg and Q m 
as shown in Figure [3] Fitting to the maximum likelihood region 
we find as(§^) ' 59 = 0.84 ± 0.05, where the quoted error is 
la for a hard prior of f2 m = 0.24. This result assumes a flat 
ACDM cosmology a nd ad opts the non-linear matter power spec- 
trum of ISmith et al.l d2003h . The redshift distribution is estimated 
from the high confidence CFHTLS-Deep photometric redshifts us- 
ing the standard n(z) model given by Eq.JlOt. Marginalisation was 
performed over h 6 [0.64, 0.80] with flat priors, n(z) with Gaus- 
sian priors as described in §6.1, and a calibration bias of the shear 
signal with flat priors as discussed in §4. The corresponding con- 
straints from each survey are presented in Figure |4]and tabulated in 
Table [4] Results in Table [4] are presente d for two methods of cal- 
cul ating the non-linear power spectrum; IPeacock & Doddsl d 19961) 
an d lSmithetal.ld2003h . We find a difference o f approximately 3% 
in the best-fit erg values, where results using ISmith et al] < t2003h 
are consistently th e smaller of the two. The more recent results of 
ISmith et ak 1 2003 ) are mor e accur ate, and should be preferred over 
those of lPeacock & Doddsl Jl996» . 

The contribution to the error budget due to the conserva- 
tive addition of the B-modes to the covariance matrix is small, 
amounting to at most 0.01 in the la error bar for a hard prior of 
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Figure 5. Values for erg when f2 m is taken to be 0.24, filled circles (solid) 
give our results with 1<t error bars, open circles (dashed) show the results 
from previous analyses (Table [TJ. Our results are given for the high con- 
fidence photometric redshift catalogue, using the functional form for n(z) 
given by Eq.(T(j), and the Smith et al. 1 2003) prescription for th e non-linear 
power spectrum. The literature values use the lPeacock & Doddsl fl9961) pre- 
scription for non-line ar power, an d are e xpected to be ~ 3% higher than 
would be the case for lSmith et alj J2003I) . The forward slashed hashed re- 
gion (enclosed by solid lines) shows the la range allowed by our combined 
result, the back slashed hashed region (enclosed by dashed lines) shows the 
1<t range given by the WMAP 3 year results. 
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Table 4. Joint constraints on the amplitude of the matter power spectrum <rg, and the matter energy density f2 m , parametrised as crg( g-jj)" = /3. The errors 
are la limits, calculated for a hard prior of f2 m = 0.24. 



Eq.{T0} Eq.OTJ Eq.(l0) Eq.(TTJ 

0a/3a/3a/3a 

CFHTLS-Wide 0.84 ± 0.06 0.55 0.81 ± 0.07 0.54 0.86 ± 0.06 0.56 0.84 ± 0.06 0.55 

GaBoDS 0.93 ±0.08 0.60 0.89 ± 0.09 0.59 1.01 ±0.09 0.66 0.98 ±0.10 0.63 

Smith et al. (20031 RCS 0.75 ± 0.07 0.55 0.73 ± 0.08 0.55 0.78 ±0.07 0.57 0.76 ±0.07 0.56 

VIRMOS-DESCART 0.99 ± 0.08 0.59 0.95 ± 0.07 0.58 1.02 ±0.07 0.60 1.00 ±0.08 0.59 

Combined 0.80 ± 0.05 0.57 0.77 ±0.05 0.56 0.84 ± 0.05 0.59 0.82 ± 0.06 0.58 



CFHTLS-Wide 
GaBoDS 
Peacock & Dodds (1996) RCS 

VIRMOS-DESCART 
Combined 



0.87 ±0.07 0.57 

0.96 ±0.09 0.62 

0.77 ±0.07 0.57 

1.02 ±0.07 0.60 

0.82 ±0.05 0.58 



0.83 ±0.06 0.56 

0.93 ± 0.08 0.60 

0.74 ±0.07 0.56 

0.98 ± 0.07 0.59 

0.79 ±0.05 0.57 



0.89 ±0.06 0.58 

1.04 ±0.09 0.66 
0.81 ±0.08 0.58 

1.05 ±0.07 0.62 
0.86 ±0.06 0.60 



0.87 ±0.06 0.57 

1.01 ±0.09 0.64 

0.79 ±0.08 0.58 

1.03 ±0.08 0.61 

0.84 ±0.06 0.59 



0, m = 0.24. The analysis was also performed having removed all 
scales where the B-modes are not consistent with zero (see Fig- 
ure[TJ. The resulting change in the best fit cosmology is small, shift- 
ing the best-fit as for an Q m of 0.24 by at most 0.03 (for CFHTLS- 
Wide), and on average by ~ 0.01 across all four surveys. These 
changes are well within our error budget. 

We present a comparison of the measured as values with those 
previously published in Figure [5] The quoted erg values are for a 
vertical slice through as — fi m space at an Sl m of 0.24, and all 
error bars denote the la region from the joint constraint contours. 
Our error bars (filled circles) are typically smaller than those from 
the literature (open circles) mainly due to the improved estimate of 
the redshift distribution. Our updated result for each survey agrees 
with the previous analysis within the error bars, this remains true 
if the literature values are lowered by ~ 3% to account for the 
difference in methods used for the non-linear power specrum. Also 
plotted are the la limits for the combined result and WMAP 3 year 
constraints (forward and back slashed hash regions respectively), 
our result is consistent with WMAP at the la level. 

In addition to our main analysis we have investigated the im- 
pact of using four different models for the redshift distribution. 
These models are dependent on which functional form is used 
(Eq.lllOt or Eg. i ll It ) and which range is used for the photometric 
redshift sample (0.0 «S z p ^ 4.0 or 0.2 < z p ^ 1.5). Tableg]gives 
the best fit joint constraints on as and f2 m , for each survey as well 
as the combined survey result. 

Comparing the best fit as values with the average redshifts 
listed in Table [2] we find, as expected, that higher redshift models 
result in lower values for as- Attempting to quantify this relation 
in terms of mean redshift fails. Changes in mean redshift are large 
(~ 14%) between the photometric samples, compared to ~ 6% 
between the two n(z) models, however a ~ 5% change in as is 
seen for both. The median redshift is a much better gauge, chang- 
ing by ~ 4% between both photometric samples and n(z) models. 
Precision cosmology at the 1% level will necessitate roughly the 
same level of precision of the median redshift. Future cosmic shear 
surveys will require thorough knowledge of the complete redshift 
distribution of the sources, in particular to what extent a high red- 
shift tail exists. 



7 DISCUSSION AND CONCLUSION 

We have performed an analysis of the 100 square degree weak lens- 
ing survey that combines four of the largest weak lensing datasets in 
existence. Our results provide the tightest weak lensing constraints 
on the amplitude of the matter power spectrum erg and matter den- 
sity fl m and a marked improvement on accuracy compared to pre- 
vious results. Using th e non-linear predic tion of the cosmological 
power spectra given in lSmith et alj d2003h . the high confidence re- 
gion of the photometric redshift calibration sample, and Eq.JlOt to 
model the redshift distribution, we find erg (£^) ' 59 = 0.84±0.05 
for a hard prior of Q m = 0.24. 

Our analysis differs from previous weak lensing analy- 
ses in three important aspects. We correctly a ccount for non- 
Gauss ian sample variance using the method of ISemboloni et aD 
( 120071) , t hus improving upon th e purely Gaussian contribution 
jven by Schneider et al.l 12002). Using the results from STEP 



Massev et alll2007al) we correct for the shear calibration bias and 
marginalise over our uncertainty in this correction. In addition we 
use the largest d eep photometric redshift catalogue in existence 
dllbert et all l2006t) to provide accurate models for the redshift dis- 
tribution of sources; we als o account for the effects of sample vari- 
ance in these distributions jVan Waerbeke et al.L[200r3) . 

Accounting for the non-Gaussian contribution to the shear 
covariance matrix, which dominates on small scales, is very im- 
portant. The non-Gaussian contribution is about twice that of the 
Gaussian contribution alone at a scale of 10 arcminutes, this dis- 
crepancy increases to about an order of magnitude at 2 arcmin- 
utes. This increases the errors on the shear correlation function 
at small scales, leading to slightly weaker constraints on cosmol- 
ogy. Ideally we would estimate the shear covariance directly from 
the data, as is done in th e previou s analyses for both G aBoDS 
faetterscheidt et al.Ll2006h and RCS faoekstra et aill2002 <J) . How- 
ever, since we can not accomplish this for all the surveys in this 
work, due to a deficit of independent fields, we opt for a consis- 
tent approach by using the analytic treatment along with th e non- 
Gaussian calibration as described in lSemboloni et all j2007l) . 

Accurate determination of the redshift distribution of the 
sources is crucial for weak lensing cosmology since it is strongly 
degenerate with cosmological parameters. Past studies using 
small external photometric catalogues (such as the HDF) suf- 
fered from sample variance, a previously unquantified source of 
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Figure 6. Redshift distri butions for the RCS survey. The solid line shows 
the best fit n(z) from lHoekstra et all l2002al) . the dashed curve is our best 
fit n(z), the average redshifts are 0.6 and 0.78 respectively. 



error that is taken into a ccount here using the prescription of 
IVan Waerbeke et all J2006h . In most cases the revised redshift dis- 
tributions were in reasonable agreement with previous results, how- 
ever this was not the case for the RCS survey. We found that the 
previous estimation was biased toward low redshift (see Figure[6]l, 
resulting in a significantly larger estimation of a$ (0.861q'q5 for 
fi m = 0.3) than is presented here (0.69 ± 0.07). This difference 
is primarily due to a problem in the filter set conversion, between 
those used to image HDF (F300W, F450W, F606W, and F814W) 
and the Cousins Rc filter used by RCS. The less severe changes 
in as for the other surveys, shown in Figure [5] can also be largely 
understood from the updated redshift distributions. Comparing the 
median redshifts from the previous results (Table Q} to those pre- 
sented here (Table[2} there exists a clear trend associating increases 
in median redshift to decreases in estimated as (see Figure |5) and 
vice-versa. Our estimate of median redshift for CFHTLS-Wide in- 
creases by ~ 10% and a corresponding decrease in as is seen, a 
similar correspondence is seen for VIRMOS-DESCART where the 
decreased median redshift in this analysis results in a proportion- 
ately larger value of us. The change in as for GaBoDS, however, 
can not be understood by comparing median redshifts since the pre- 
vious results use a piece wise function to fit the redshift distribution, 
hence comparing median redshifts is not meaningful. 

Our revised cosmological constraints introduce tension, at the 
~ 2a level, between the results from the VIRMOS-DESCART 
and RCS survey as shown in Figure [5] In this analysis we have 
accounted for systematic errors associated with shear measure- 
ment but have neglected potential systematics arising from cor- 
relations between galaxy shape and the underlying density field. 
This is valid in the case of intrinsic galaxy alignments which are 
expected to contribute less than a percent of th e cosmic shear 
signal for the deep surveys used in this analysis dHevmans et all 
2006b). What is currently uncertain however is the level of 
systematic error that arises from shear-ellipticity correlations 
dHirata & Seliakll2004lHevmans et afll2006bl.lMandelbaum et aT 



2006; Irlirata et aUl2007l) which could reduce the amplitude of the 
measu red shear correlation function by ~ 10% dHevmans et all 
2006b). It h as been found from the analysis of the Sloa n Digital 
Sky Survey dlvlandelbaum et alll2006l ; lHirata et alll2007l) . that dif- 
ferent morphological galaxy types contribute differently to this ef- 
fect. It is therefore possible that the different R-band imaging of the 
RCS and the slightly lower median survey redshift make it more su- 



ceptible to this type of systematic error. Without complete redshift 
information for each survey, however, it is not possible to test this 
hypothesis or to correct for this potential source of error. In the 
future deep multi-colour data will permit further investigation and 
correction for this potential source of systematic error. 

In addition to using the best available photometric redshifts, 
we marginalise over the redshift distribution by selecting parame- 
ter triplets from their full 3D probability distribution, instead of fix- 
ing two parameters and varying the third as is often done. Thus the 
marginalisation is representative of the full range of n(z) shapes, 
which are difficult to probe by varying one parameter due to degen- 
eracies. 

We have conducted our analysis using two different func- 
tional forms for the redshift distribution, the standard form given 
in Eq.dTUt and a new form given by Eg. ill It. The new function is 
motivated by the presence of a high-z tail when the entire photo- 
metric catalogue (0.0 ^ z p ^ 4.0) is used to estimate the redshift 
distribution, Eq.dlOt does a poor job of fitting this distribution (Fig- 
ure |2j- However when we restrict the photometric catalogue to the 
high confidence region (0.2 ^ z p ^ 1.5) both functions fit well, 
and the tendency for Eq.il It to exhibit a tail towards high-z in- 
creases the median redshift resulting in a slightly lower estimate of 
as- Though the tail of the distribution has only a small fraction of 
the total galaxies, they may have a significant lensing signal owing 
to their large redshifts. The influences of the different redshift dis- 
tributions on cosmology is consistent within our la errors but as 
survey sizes grow and statistical noise decreases, such differences 
will become significant. As the CFHTLS-Deep will be the largest 
deep photometric redshift catalogue for some years to come, this 
posses a serious challenge to future surveys attempting to do preci- 
sion cosmology. To assess the extent of any high redshift tail, future 
surveys should strive to include photometric bands in the near infra- 
red, allowing for accurate redshift estimations beyond z = 1.5. 

For fi m = 0.24, the combined results using Eq.lllOt are 
0.80 ± 0.05 and 0.84 ± 0.05 for the unrestricted and high con- 
fidence photometric redshifts respectively, using Eq.il It we find 
0.77± 0.05 and 0. 82 ±0.06 (these res ults use the non-linear power 
spectrum given by Smith et al. (2003)). For com pleteness we pro- 
vide re sults for both the lSmith etall d2003l) and lPeacock & D odds 
( 1996) non-li near power spectra, the resulting as values differ by 
~ 3% of the lSmithetai1d2003l) v alue w hich is consistently the 
smaller of the two. The Smi th et al. | d2003h study is known to pro- 
vi de a more accurate estima tion of the non-linear power than that 
of iPeacock & Dodds! d 1996b. for this reason we prefer the results 
obtained using the lSmith et all ((2003) model. However, given the 
magnitude of variations resulting from different redshift distribu- 
tions, this difference is not an important issue in the current work. 
Accurately determining the non-linear matter power spectrum is 
another challenge for future lensing surveys intent on precision 
cosmology. Alternatively, surveys focusing only on the large scale 
measurement of the shear (Fu et al., in prep) are able to avoid com- 
plications arising from the estimation of non-linear power on small 
scales. 

Surveys of varying depth provide different joint constraints 
in the fi m — erg plane, thus combining their likelihoods produces 
some degeneracy breaking. Taking the preferred prescription for 
the non-linear power spectrum, and using Eq.lllOt along with the 
high confidence photometric redshifts to estimate the redshift dis- 
tribution, we find an upper limit of fl m <~ 0.4 and a lower limit 
of as >~ 0.6 both at the ler level (see Figure[3]l. 

Weak lensing by large scale structure is an excellent means 
of constraining cosmology. The unambiguous interpretation of the 
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shear signal allows for a direct measure of the dark matter power 
spectrum, allowing for a unique and powerful means of constrain- 
ing cosmology. Our analysis has shown that accurately describing 
the redshift distribution of the sources is vital to future surveys in- 
tent on precision cosmology. Including near-IR bands to photomet- 
ric redshift estimates will be a crucial step in achieving this goal, 
allowing for reliable redshift estimates at z > 1.5 and assessing the 
extent of a high-z tail. 
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Table Al. E and B modes of the shear correlation function for the CFHTLS- 
Wide survey, the error (<5£) is statistical only, and given as the standard 
deviation. 



8 (arcmin) 


He 


Hb 


SH 


0.74 


1.53060e-04 


-4.16557e-05 


2.31000e-05 


1.44 


1.00161e-04 


1.79443e-05 


1.22000e-05 


2.37 


5.82357e-05 


-1.60557e-05 


9.64000e-06 


3.54 


3.11136e-05 


7.14429e-06 


6.66000e-06 


5.41 


2.89396e-05 


-2.05571e-06 


4.29000e-06 


7.28 


3.08030e-05 


-3.47571e-06 


4.76000e-06 


8.69 


1.41053e-05 


6.14429e-06 


4.43000e-06 


11.72 


1.66208e-05 


-5.15714e-07 


2.16000e-06 


18.74 


1.50680e-05 


4.14286e-07 


1.25000e-06 


28.09 


1.07719e-05 


-7.957 15e-07 


1.03000e-06 


37.44 


8.72227e-06 


4.14286e-07 


9.10000e-07 


49.12 


8.11150e-O6 


-2.18571e-06 


6.64000e-07 


62.92 


8.97071e-06 


-3.01571e-06 


6.2500()e-()7 



Table A2. E and B modes of the shear correlation function for GaBoDS, 
the error (<5£) is statistical only, and given as the standard deviation. 



8 (arcmin) 


He 


Hb 






0.36 


2.80637e-04 


1.38754e-04 


9.54000e- 


•05 


0.50 


1.72930e-04 


1.24754e-04 


6.92000e- 


•05 


0.70 


1.15363e-04 


6.74540e-05 


4.99000e- 


■05 


0.98 


1.20006e-04 


-1.59760e-05 


3.81O00e- 


■05 


1.36 


8.75079e-05 


2.76540e-05 


2.8()000e- 


■05 


1.90 


7.4323 le-05 


8.35400e-06 


2.15000e- 


■05 


2.65 


9.55859e-05 


-2.30860e-05 


1.71000e- 


■05 


3.69 


6.48524e-05 


3.05400e-06 


1.40000e 


■05 


5.16 


4.30325e-05 


5.25400e-06 


1.19000e- 


■05 


7.19 


2.73408e-05 


-7.24600e-06 


1.03000& 


■05 


10.03 


9.97772e-06 


-7.46000e-07 


8.65000e- 


■06 


13.99 


1.22154e-05 


-6.54600e-06 


6.55000e- 


■06 
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APPENDIX A: SHEAR CORRELATION FUNCTION AND 
COVARIANCE MTRIX FOR THE SURVEYS 

We present both the measured shear correlation function and the 
covariance matrix for each survey. The shear correlation function, 
given in Tables IaTI lA2l |A3l and lA4l for the CFHTLS-Wide, Ga- 
BoDS, RCS and VIRMOS-DESCART surveys respectively, has 
been calibrated on large scales where (Mf p (A8))± is consistent 
with zero as described in §4. 

Tables |A31 lA6l lA7l and lA8l for the CFHTLS-Wide, GaBoDS, 
RCS and VIRMOS-DESCART surveys respectively tabulate the 
correlation coefficient matrix: 



J ~ <a i/2 ^> i/2 ' 

where C is the covariance matrix for each survey, as described in 
§6. We also tabulate (£f ) so that the covariance matrix may be cal- 
culated from the correlation coefficient matrix. Note that (£, 2 ) is the 
variance of the i th scale, equivalent to Cu. 



Table A3. E and B modes of the shear correlation function for RCS, the 
error (<5£) is statistical only, and given as the standard deviation. 



8 (arcmin) 


He 


Hb 


SH 




1.12 


9.71718e-05 


2.78320e-05 


1.99000e- 


■05 


1.75 


3.38178e-05 


1.67320e-05 


1.91000e- 


■05 


2.50 


3.74410e-05 


4.00200e-06 


1.16000e- 


■05 


3.24 


4.27205e-05 


6.63200e-06 


1.43000e- 


■05 


4.00 


1.49772e-05 


1.58320e-05 


9.21000e- 


06 


5.75 


1.52878e-05 


-8.80000e-08 


5.13000e- 


06 


8.50 


1.50807e-05 


-3.89600e-06 


3.95000e- 


■06 


12.25 


6.89234e-06 


2.15200e-06 


2.85000e- 


■06 


16.49 


4.81573e-06 


1.83200e-06 


2.65000e- 


■06 


20.74 


4.6915 le-06 


1.82200e-06 


2.30000e- 


■06 


26.49 


5.29607e-06 


3.59200e-06 


1.75000e- 


■06 


34.99 


5.33747e-06 


1.11200e-06 


1.36000e- 


■06 


46.49 


6.91304e-06 


-4.88000e-07 


1.08000e- 


■06 


56.51 


6.69565e-06 


-4.88000e-07 


1.34000e- 


■06 



Table A4. E and B modes of the shear correlation function for the 
VIRMOS-DESCART survey, the error (<5£) is statistical only, and given as 
the standard deviation. 



8 (arcmin) 


He 


Hb 


SH 


0.73 


2.26695e-04 


3.91133e-05 


3.98O00e-05 


1.05 


1.69759e-04 


-1.59667e-06 


3.90000e-05 


1.51 


1.30421e-04 


5.31333e-06 


2.20000e-05 


2.29 


9.60525e-05 


2.82133e-05 


1.57000e-05 


3.20 


6.603 18e-05 


-1.38767e-05 


1.35000e-05 


4.11 


5.56798e-05 


6.41333e-06 


1.22000e-05 


5.02 


5.29883e-05 


3.71333e-06 


1.12000e-05 


6.39 


3.50794e-05 


-1.26667e-06 


7.46000e-06 


8.21 


1.94893e-05 


2.11333e-06 


6.75000e-06 


11.41 


2.721 19e-05 


-8.46667e-07 


3.99000e-06 


15.97 


1.98620e-05 


4.81333e-06 


3.50000e-06 


20.53 


1.67357e-05 


5.41333e-06 


3.23000e-06 


25.25 


1.05659e-05 


4.51333e-06 


2.95000e-06 


30.23 


9.77399e-06 


1.31333e-06 


2.73000e-06 


35.10 


1.06073e-05 


-1.05667e-06 


2.71000e-06 


39.78 


1.1 777 le-05 


-2.08667e-06 


2.61000e-06 
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Table A5. The correlation coefficient matrix (r) for the CFHTLS-Wide survey, the scales 0; correspond to those given in Table lAll The column (£ 2 ) is given 
in units of 10~ 10 , and can be used to reconstruct the covariance matrix C. 







01 


02 


03 


04 


05 


06 


07 


08 


09 


010 


0u 


012 


013 


01 


29.9246 


1.0000 


0.2595 


0.1790 


0.1975 


0.1900 


0.1223 


0.0841 


0.1025 


0.0707 


0.0587 


0.0549 


0.0393 


0.0266 


e-z 


6.5960 


0.2595 


1.0000 


0.2175 


0.2387 


0.2330 


0.1519 


0.1062 


0.1527 


0.1431 


0.1251 


0.1169 


0.0837 


0.0567 


03 


4.2569 


0.1790 


0.2175 


1.0000 


0.2182 


0.2192 


0.1585 


0.1243 


0.1904 


0.1784 


0.1560 


0.1457 


0.1044 


0.0707 


04 


1.3550 


0.1975 


0.2387 


0.2182 


1.0000 


0.3734 


0.2861 


0.2228 


0.3405 


0.3180 


0.2778 


0.2594 


0.1855 


0.1254 


6r, 


0.5285 


0.1900 


0.2330 


0.2192 


0.3734 


1.0000 


0.4829 


0.3683 


0.5554 


0.5091 


0.4370 


0.4120 


0.2970 


0.2006 


06 


0.6121 


0.1223 


0.1519 


0.1585 


0.2861 


0.4829 


1.0000 


0.3505 


0.5277 


0.4732 


0.3984 


0.3790 


0.2757 


0.1865 


07 


0.8166 


0.0841 


0.1062 


0.1243 


0.2228 


0.3683 


0.3505 


1.0000 


0.4683 


0.4119 


0.3464 


0.3293 


0.2393 


0.1615 


08 


0.2533 


0.1025 


0.1527 


0.1904 


0.3405 


0.5554 


0.5277 


0.4683 


1.0000 


0.7518 


0.6322 


0.5951 


0.4289 


0.2901 


00 


0.1796 


0.0707 


0.1431 


0.1784 


0.3180 


0.5091 


0.4732 


0.4119 


0.7518 


1.0000 


0.7832 


0.7144 


0.5112 


0.3522 


m 


0.1447 


0.0587 


0.1251 


0.1560 


0.2778 


0.4370 


0.3984 


0.3464 


0.6322 


0.7832 


1.0000 


0.8248 


0.5722 


0.3883 


011 


0.1074 


0.0549 


0.1169 


0.1457 


0.2594 


0.4120 


0.3790 


0.3293 


0.5951 


0.7144 


0.8248 


1.0000 


0.6889 


0.4453 


012 


0.1274 


0.0393 


0.0837 


0.1044 


0.1855 


0.2970 


0.2757 


0.2393 


0.4289 


0.5112 


0.5722 


0.6889 


1.0000 


0.4240 


013 


0.1615 


0.0266 


0.0567 


0.0707 


0.1254 


0.2006 


0.1865 


0.1615 


0.2901 


0.3522 


0.3883 


0.4453 


0.4240 


1.0000 



Table A6. The correlation coefficient matrix (r) for GaBoDS, the scales 9i correspond to those given in Table lA2l The column (£ 2 ) is given in units of 10 
and can be used to reconstruct the covariance matrix C. 





(O 


0i 


02 


03 


04 


05 


06 


07 


08 


09 


010 


0n 


012 


01 


323.7553 


1.0000 


0.1199 


0.1298 


0.1737 


0.1383 


0.1627 


0.1028 


0.1404 


0.1148 


0.0890 


0.0868 


0.0590 


02 


223.3641 


0.1199 


1.0000 


0.1214 


0.1566 


0.1249 


0.1490 


0.0958 


0.1336 


0.1117 


0.0892 


0.0896 


0.0604 


03 


80.2861 


0.1298 


0.1214 


1.0000 


0.2051 


0.1608 


0.1937 


0.1269 


0.1810 


0.1552 


0.1284 


0.1351 


0.0881 


e± 


22.1309 


0.1737 


0.1566 


0.2051 


1.0000 


0.2487 


0.2988 


0.1991 


0.2902 


0.2548 


0.2150 


0.2286 


0.1514 


0-, 


18.2583 


0.1383 


0.1249 


0.1608 


0.2487 


1.0000 


0.2817 


0.1891 


0.2804 


0.2536 


0.2193 


0.2327 


0.1654 


00 


6.9091 


0.1627 


0.1490 


0.1937 


0.2988 


0.2817 


1.0000 


0.2815 


0.4205 


0.3819 


0.3349 


0.3720 


0.2681 


07 


9.2132 


0.1028 


0.0958 


0.1269 


0.1991 


0.1891 


0.2815 


1.0000 


0.3416 


0.3155 


0.2839 


0.3190 


0.2314 


9s 


2.6680 


0.1404 


0.1336 


0.1810 


0.2902 


0.2804 


0.4205 


0.3416 


1.0000 


0.5807 


0.5259 


0.5879 


0.4247 


00 


2.1958 


0.1148 


0.1117 


0.1552 


0.2548 


0.2536 


0.3819 


0.3155 


0.5807 


1.0000 


0.5801 


0.6401 


0.4606 


010 


2.0117 


0.0890 


0.0892 


0.1284 


0.2150 


0.2193 


0.3349 


0.2839 


0.5259 


0.5801 


1.0000 


0.6561 


0.4734 


en 


1.1340 


0.0868 


0.0896 


0.1351 


0.2286 


0.2327 


0.3720 


0.3190 


0.5879 


0.6401 


0.6561 


1.0000 


0.6095 


012 


1.1960 


0.0590 


0.0604 


0.0881 


0.1514 


0.1654 


0.2681 


0.2314 


0.4247 


0.4606 


0.4734 


0.6095 


1.0000 



Table A7. The correlation coefficient matrix (r) for RCS, the scales 8i correspond to those given in Table lA3l The column {^) is given in units of 10 
and can be used to reconstruct the covariance matrix C. 





(O 


01 


02 


03 


04 


05 


06 


07 


08 


09 


010 


0ii 


012 


013 


014 


01 


12.9511 


1.0000 


0.0957 


0.1251 


0.0775 


0.0543 


0.1070 


0.0731 


0.0695 


0.0618 


0.0578 


0.0416 


0.0535 


0.0484 


0.0389 


02 


7.0289 


0.0957 


1.0000 


0.1354 


0.0822 


0.0574 


0.1155 


0.0832 


0.0885 


0.0832 


0.0784 


0.0564 


0.0725 


0.0656 


0.0526 


03 


1.8431 


0.1251 


0.1354 


1.0000 


0.1445 


0.0996 


0.2076 


0.1575 


0.1721 


0.1626 


0.1529 


0.1100 


0.1414 


0.1278 


0.1024 


04 


2.7380 


0.0775 


0.0822 


0.1445 


1.0000 


0.0788 


0.1696 


0.1290 


0.1407 


0.1333 


0.1254 


0.0902 


0.1159 


0.1048 


0.0840 


05 


3.5622 


0.0543 


0.0574 


0.0996 


0.0788 


1.0000 


0.1535 


0.1149 


0.1229 


0.1161 


0.1099 


0.0791 


0.1016 


0.0919 


0.0736 


00 


0.4239 


0.1070 
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Table A8. The correlation coefficient matrix (r) for the VIRMOS-DESCART survey, the scales 0; correspond to those given in Table lA4l The column (£ 2 ) is 
given in units of 10 -10 , and can be used to reconstruct the covariance matrix C. 
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